Mott insulating state in a quarter-filled two-orbital 
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Abstract. We investigate the ground-state properties of the one-dimensional two-band 
Hubbard model with different bandwidths. The density-matrix renormalization group method is 
applied to calculate the averaged electron occupancies n as a function of the chemical potential /j,. 
Both at quarter and half fillings, "charge plateaux" appear in the n-/j plot, where dn/dn diverges 
and the Mott insulating states are realized. To see how the orbital polarization in the one-quarter 
charge plateau develops, we apply the second-order perturbation theory from the strong-coupling 
limit at quarter filling. The resultant Kugel-Khomskii spin-orbital model includes a magnetic 
field coupled to orbital pseudo-spins. This field originates from the discrepancy between the 
two bandwidths and leads to a finite orbital pseudo-spin magnetization. 



' "Orbital degrees of freedom" has been one of the important keywords to understand the 

low temperature properties of strongly correlated electron systems. Among them, recently, the 



orbital-selective Mott transition (OSMT) is proposed to explain the exotic heavy metallic state in 
Ca2_xSra;Ru04[Tl[2]. To illustrate this scenario, the two-orbital model with different bandwidths 
OO ! has been extensively investigated by dynamical mean-field theory |3l !!]• They have clarified the 

■ realization conditions of the OSMT by changing various parameters, such as Coulomb interaction 

parameters, band structures, band filling, and others. In this paper, we also study the ground 
state properties of the two-band Hubbard model with different bandwidths, but in one dimension 
(ID). Our main concern here is whether the OSMT survives in ID, where the effect of quantum 
5^ i fluctuations is most severely enhanced. 

First, let us begin with the description of the ID two-orbital Hubbard model[5j defined by; 

ijacr ioL iaa' 

-jY.4laCil<r'42a'Ci2<T - jJ2{4l^4nCi2]Ci2l + h.C.^ , (1) 
iaa' i 

where c^^,^ creates an electron with spin a (=1 or [) and orbital a (= 1 or 2) at the j-th site. 
The electron transfers of strength ta are possible between the same type of neighboring orbitals 
and /X denotes the chemical potential. As for the intra-site Coulomb interactions, we assume 



the rotational invariance for simplicity, and thus the Coulomb parameters, U , U', and J, always 
hold the following relation; U = U' + 2J. As long as the intra-atomic Coulomb interactions 
are concerned, physically relevant parameters would satisfy U > U' > J. However, we do not 
restrict our calculations to this parameter region, since in general the above Hamiltonian can 
be also viewed as a variation of coupled Hubbard chains. 

In this paper, we have used the infinite-size density-matrix renormalization group (DMRG) 
algorithm[ni [7] modified by adapting the recursion relation to the wave function, namely, 
the wave function is so upgraded that it should gradually approach the ground-state wave 
function when the system size is extended. This technique, called the product-wave-function 
renormalization-group method [SI E]) enhances the accuracy of calculations and overcomes 
metastability problems during the renormalizing process if the ground state to target is gapless. 
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Figure 1. (Color online) The averaged electron densities of orbital 1 (ni, denoted by open 
circles), orbital 2 (n2, solid triangles) and their sum (ntot; open squares) as a function of chemical 
potential /x. Coulomb interactions are set to be C/ = 10, U' = 5 and J = 2.5. 

The different bandwidths are implemented in Eq.([T]) by taking ti = 1 and t2 = 0.5. Since we 
are mainly interested in strongly correlated regime, U is always fixed at 10 throughout this paper 
and U' as well as J are varied while keeping U = U' + 2 J. By means of the DMRG technique, 
we have calculated the electron density of a-orbital defined by Ua = J2ia{^iacT) / ^ , where N 
is the total number of sites. Figure 1 shows the ;U-dependence of ni, 71-2, and ntot = + 7^2 
for U' = 5 and J = 2.5. Let us take a closer look at Fig. [1] from the low chemical potential 
region. When /u exceeds /id = —2ti, the bottom of the orbital-1 energy band, ni increases 
sharply while n2 stays zero until /i reaches iic2- At fic2, ni and 712, respectively, drops down 
and jumps up suddenly. Then between fj,c2 and /i^s, ni (712) shows a gradual (steep) monotonic 
increase and, finally, the first plateaux for both ni and n2 appear simultaneously between fj,qs 
and /iqe- At these plateaux, the values of ni and n2 are irrational but ntot definitely equals 
to unity corresponding to one-quarter filling. Therefore, this Mott insulating state itself is not 
controversial except for that ni is not equal to n2- We will come back to this point later. With 
further increase of n above iJ,qe in Fig. [H we find the next plateau at n2 = 1 between ij,hs2 and 
^hsi, which means the charge gap opens for the orbital-2 band. On the other hand, ni increases 
continuously toward unity in this region and thus the orbital-1 band should be metallic. This 
implies that the filling-control OSMT takes place at fihs2- Then, the orbital-selective Mott state 
with one orbital localized and the other itinerant is realized just below one-half filling in the 
present system. For /i > fihsi, the ni = 1 plateau emerges in addition to the n2 = 1 plateau, 
which correspond to the normal half-filled Mott insulating state. 
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Figure 2. (Color online) The averaged electron densities, ni and n2 (upper panel), and the 
nearest neighbor spin-spin and orbital-orbital correlation functions (lower panel) as a function 
of U', where ntot = 1, U = 10 and U = U' + 2 J holds. 



Hereafter, we will investigate the quarter-filled Mott state in more detail. For this 
purpose, the DMRG calculations are performed within the S^^^ = J2i Sf = subspace at 
quarter filling. In Fig. [2] shown are the electron densities (ni and n2, upper panel) and 
correlation functions {{TfT-i^i) and {Si ■ 5'j+i), lower panel) as a function of U' for U = W 
and J = (U — U')/2. Considering the atomic limit, there are three possible ground states 
characterized by; i) spin-triplet and orbital-singlet states like c|]^i-c|2-|-|0) (averaged site energy 
e = ([/' — J)/2 = (3[/' — U)/4), ii) singly occupied states like c|q,^|0) (null site energy), and iii) 
spin-singlet and orbital-triplet states hke (c|;^|^c|-^j^ -|- c|2T42i)l'-*) = + = {^U — U')/A). 
Therefore, two quantum critical points, C/^i = U/3 and = 3C/, separate the above three 
phases when ^1=^2 = 0. The singularity around U' = 3.5 in Fig. [2] seems to correspond 
to U^i and, in fact, the vanishing orbital polarization below U' < 3.5 is consistent with 
the formation of a local orbital singlet. With further increase of U' above U^i, an orbital 
magnetization: (Tf) = (ni — n2)/2 develops accompanied by the antiferromagnetic spin-spin 
and the ferromagnetic orbital-orbital correlations. These results show that, at the one-quarter 
plateau, the magnitude of charge disproportionation between rii and n2 depends on Coulomb 
interaction parameters and only the total electron density is preserved. For U' > U = 10, 
{Tf) is saturated and all electrons reside in the orbital 1. In such a strong-coupling regime, 
the system should be identical to the isotropic Heisenberg spin chain, which is confirmed by 
{Si ■ Si+i) = 1/4 - In 2 ~ -0.443 as well as {TfTf^^) = 1/4 in Fig. [2l 

In order to examine how the orbital polarization at quarter filling is developed with increasing 
U', we have derived the effective spin-orbital Hamiltonian of Kugel-Khomskii tvpe|10j. Starting 
from the atomic limit and within U^^ ^ U' < the effective Hamiltonian is given by; 

Won = T.I [wg, - aSJJ.,(3;' + r4i)] with; 
<' = -(l7^ + ^)(i-^.-^.)K^ + *=^)(i + ^.%')} 



-JPTl {\-^'' + (i - ^''^O + + ^i"^*) } 

where the spin-1/2 {S) and orbital pseudospin-1/2 (T) operators are, respectively, defined 

by Sj=J2aaa' i^'^]aa'^uu' Cjaa') /2 and fj=Y,aaa' ^^ja' o) /2 with the USe of the Pauh 

^ {S) 

matrices r's. is an effective crystal field coupled to orbital pseudospins, which originates 

from the discrepancy between the two bandwidths, that is, ti 7^ t2. Since the {U + J)~^ term 

in A-'^\ which is proportional to (3C/ — U')~^, increases with an increase of U' {U' < 

a ferromagnetic (T-^T-f,_j^) is expected to be developed accordingly. Therefore, the different 

bandwidth in the present model is indispensable for producing the orbital polarization. On 

i s) 

the other hand, when ti = t2 and U = U', A-j vanishes and TCcS is equivalent to the SU(4) 
spin-orbital model [TTlll2j . In this model, the ground state is an SU(4) singlet ground state with 
(Tf ) = corresponding to rii = n2- 

To conclude, we have investigated the ground-state properties of the ID two-orbital Hubbard 
model with different transfer integrals by the DMRG method. Two plateaux appear in the 
total electron density as a function of the chemical potential, which suggests the existence of 
Mott insulating states. The first insulating state ranging in /x^s < n < Hqe corresponds to the 
quarter-filled Mott insulator and the charge disproportionation between two orbitals exists due 
to the different bandwidth of the two bands. On the other hand, in the second Mott insulating 
state for fihsi < /x, the two bands are both half filled and thus the total electron density is at 
half filling. It is also found that the OSMT takes place just below half filling, where the charge 
gap exists in the narrower band and the wider band remains metallic. The detailed analysis 
around the charge one-half plateau will be published elsewhere. 
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